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Abstract 

We study the stability of supersymmetry breaking metastable vacua and su- 
persymmetric vacua in the presence of sohtons. The metastable vacua of su- 
persymmetric QCD and those found elsewhere such as in models based on the 
SU(5) grand unified group support the existence of topological solitons. The 
vacua containing such topological defects can become unstable against decay 
into lower energy configurations. We show for a specific model that a finite 
region of the available parameter space of couplings becomes disallowed due to 
the presence of monopoles. In a manner similar to previous studies based on 
cosmic strings, it is shown that soliton solutions arising in supersymmetric the- 
ories can put constraints on the range of allowed values of the couplings arising 
in the theories. Implications for cosmology are discussed. 

Keywords: topological soliton, supersymmetry breaking, metastable vacua, 
cosmology 
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1. Introduction 

Dynamical supersymmetry breaking in metastable vacua is an effective means 
of breaking supersymmetry which can be accommodated in various classes of 
models such as A'" = 1 supersymmetric SU{Nc) QCD with Nf massive funda- 
mental flavors [l|, Q ■ Models based on this method can have several supersym- 
metry breaking false vacua which have a finite lifetime based on the quantum 
tunneling rate to the true vacuum |3|]. Metastable vacua which break super- 
symmetry also occur in many other models of supersymmetry breaking and 
mediation 0, S H 0| ■ Issues relating to the phenomenological iinplementation 
of these ideas have been discussed in [1, [iJ3,E [iLM El, llJ, H El E El, E3| 



and cosmological applications appear in [18|, |19|, |20| ■ A necessary condition for 
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the viability of such models is that the lifetime of the metastable vacua is much 
larger than the age of the universe. 



In a recent paper |21| . we had emphasized that it is not sufficient to study 
the stability of supersymmetry breaking vacua in terms of their translation 
invariant form alone. In general, topological defects can form in the evolution 
of the early universe [2^ . For example, the models of supersymmetry breaking 
described in [l|, Q break SUSY in metastable vacua which can contain cosmic 
strings . Such defects have important consequences with regards to stability 
issues. The core of a cosmic string or monopole can give rise to a "seeding" 
of the true vacuum and render it unstable against decay into a lower energy 
configuration [24, 

M,^- This is true even if the translation invariant vacuum 
is sufficiently stable against quantum mechanical tunneling to the true vacuum. 
Such a process has been studied in the context of phase transitions in Grand 
Unified Theories in Instabilities of domain walls in theories with compact 
extra dimensions are discussed in [isl ]. 

In 21| we dealt with the messenger sector of a model in which supersymme- 
try breaking occurs according to the gauge mediated supersymmetry breaking 
scenario. The seeding resulted from the presence of cosmic string solutions. 
As an extension of this work, we study the seeding mechanism in the presence 
of monopoles, and this is done for the case of direct supersymmetry breaking 
rather than having supesymmetry breaking being communicated from a hid- 
den sector. In particular, we have chosen to study a non-abelian model in 
which both the gauge and supersymmetry breaking occur simultaneously in a 
O'Raifeartaigh type model through vacuum expectation values of appropriate 
Higgs scalars. After SU{5) breaking, there are no supersymmetric vacua in this 
model. There are however multiple supersymmetry breaking vacua, some of 
which are metastable. It turns out that some of these metastable and spatially 
homogeneous vacuum states become disallowed in the presence of monopoles for 
a large range of couplings occurring in the model. In addition to the instability 
which can arise in supersymmetry breaking vacua, we show in the present work 
that this effect can also take place for two supersymmetric vacua. It is possi- 
ble to study the stability numerically and also semi-analytically. Both of these 
approaches will be dealt with in this paper. 

The classical instabilities of solitons discussed here can have a significant 
impact on the evolution of the early universe. In the model we deal with, there 
is a metastable supersymmetry breaking vacuum which is undesirable from the 
point of view of phenomenology. It is possible for some causally connected re- 
gions of the universe to get trapped in this false vacuum during the early stages 
of its evolution. If no monopoles are present, the transition to the true vacuum 
can only take place via quantum tunneling. This can result in a potentially inho- 
mogeneous universe with the first-order phase transition never being completed 
"2^, 3^ . As we will show, the presence of monopoles in such metastable vacua 



can result in classical instabilities of such configurations and hence a graceful 
exit from a potentially inhomogeneous expansion. 

The next section [2] discusses the classical instabilities of solitons in general. 
The following section [3] then briefly describes the model and its vacua, while 
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section 0] deals with the monopole ansatz and equations of motion. Sections [5] 
and [6] discuss the stabihty of the vacua containing monopoles from a numerical 
and semi-analytic approach. Implications for cosmology are presented in section 
[7] followed by concluding remarks in section [51 

2. Classical Instabilities of Monopoles and Vortices 

We begin with brief overview of how vortex and monopole solutions can 
become unstable classically under certain conditions. The dissociation of SU (5) 
monopoles which will be discussed in this paper is similar to the monopole 
dissociation studied in [2^ . A general study of this phenomenon which applies 
for both vortices and monopoles is discussed in (26| . 
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Figure 1; An example of the potential to be studied in this paper. The value of t is given by 
the difference in energy densities of the false and true vacuum. The term ^ is a mass scale 
which will appear in the model discussed here. 

Let us consider the scalar potential for a field CTi as shown in figure [T] This 
potential is an example of the type of potential we will encounter in this study. 
It has a global minimum of nearly vanishing energy for which the VEV of ui is 
close to zero. There is also a false vacuum in which the value of ui is larger than 
that of the true vacuum. The term e denotes the difference in energy densities 
of the false and true vacuum and it is positive in this case. 

A topological soliton solution such as a cosmic string or magnetic monopole 
has a vanishing field strength within its core due to continuity requirements. 
Consequently, the field strength must increase from zero and approach the vac- 
uum value asymptotically. Now if a monopole is present in the false vacuum of 
figure [Jl the core of this monopole contains a region in which the field strength 
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corresponds to the true vacuum. Following the discussion of [26|, we assume 
that this region is a spherical bubble of radius Rb which has a thin boundary. 
This is referred to as the "thin-wall" approximation. Under this approximation, 
the monopole consists of three regions. The inner-most region is a spherical re- 
gion which corresponds to the true vacuum. This is followed by a "thin-wall" in 
which the field strength grows rapidly to value corresponding to the metastable 
vacuum. Outside this wall is a region of false vacuum in which the energy 
density is higher than that of the true vacuum. 

A question which naturally arises is whether such a bubble of true vac- 
uum can expand and occupy the region of false vacuum outside the core of the 
monopole. If the bubble expands from a radius R to R + dR, the change in 
energy to first order in dR is —AireR^dR. The negative sign signifies that en- 
ergy can be lost through radial expansion. Thus, the energy of the bubble is 
proportional to —eR^. However, an increase in radius of the bubble leads to 
an increase in surface area of the surrounding wall which has a positive energy 
density a. Hence, there is a term in the total bubble energy proportional to 
aR? . There is also a gauge field present in a topological soliton and its energy 
varies as C /R for some constant C. Thus, the total energy of the bubble is given 

by 

E{R) = -€^R^ + A'KaR^ + (1) 
o R 

This energy as a function of bubble radius R is plotted in figure [2j The profile 
E{R) has a local minimum for small values of e. The value of R for this minimum 
is denoted Rh and this corresponds to a stable monopole solution. As e increases, 
the curvature of the local minimum becomes less positive and at a critical value 
eo, it becomes zero. If e is increased beyond cq, there is no local minimum and 
Rb tends to infinity. This indicates an unstable monopole configuration whose 
core contains a region of true vacuum which is expanding into the surrounding 
false vacuum. 

In the supersymmetry breaking model we will consider, it turns out that 
value of e is governed by a parameter M . Increasing M increases both the VEV 
of (7i and also the energy of the false vacuum (see figure [T]). Thus, 

e - (2) 

As will be shown in a later section, increasing M beyond a certain value creates 
an instability in the monopole configuration which settles into the metastable 
vacuum. 

We will be studying instability from two approaches. One is by simply look- 
ing for time-independent monopole solutions. If a monopole becomes classically 
unstable in the way described above, it is described by a time-dependent solution 
describing an expanding bubble of true vacuum. In this case, it is impossible 
to obtain a time-independent solution. Hence, whenever a time-independent 
solution cannot be obtained, that monopole configuration is understood to be 
unstable for those values of parameters. This will be discussed in greater detail 
in section [01 We also follow a semi-analytic approach to the study of stability 
and this will be discussed in section [521 
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RADIUS OF BUBBLE -- 



Figure 2: The energy of a bubble of radius R. For small enough e, there is a minimum of E(R) 
corresponding to a stable monopole solution with finite bubble radius. If e becomes too large, 
there is no minimum and Ri, — >■ oo. This corresponds to an expansion of the core resulting in 
an unstable monopole configuration. 



3. The model and its vacua 



We study a model described in [31[ which enables metastable supersymmetry 
breaking without the use of singlet fields. This requires the use of two SU (5) 
adjoints Si and S2 whose vacuum expectation values (VEVs) are given by 
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where Ci (i = 1, 2) are the standard model singlets. The superpotential is 

= -i- az (sOMfi - VsOMXai + {7ai + 30a2)al) ai 
from which the scalar potential can be written as 



(4) 
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There are no D-Term contributions to the scalar potential since the VEVs of 
Ei are diagonal. Defining vi = (ai) and V2 — (CT2), SU{5) is broken at V2 while 
supersymmetry is broken at vi. The VEV vi can be obtained from minimizing 
but V2 is undetermined at tree order. This flat direction is lifted by nontrivial 
1-loop corrections to the scalar potential which have the form V ss m^^ (c2— ■^2)^- 
Adding a term hke this on the right hand side of (O leaves the value of fi 
unchanged while ((T2) gets fixed at V2- The units of V2 are the same as those of 
(72 and can be taken to be 100^ to simplify the forthcoming discussion. We thus 
stabilize V2 in the ensuing numerical simulations by adding ((T2 — ^2)'^ to the 
scalar potential (O and hence fixing a2 to V2 (in units of 100/i). An example of 
the scalar potential with V2 and the couplings fixed is shown in figure [T] 

3.1. Supersymmetry breaking vacua 

The model we are considering has no supersymmetric minima unless V2 = 0. 
When SU (5) breaks with V2 ^ 0, there are two supersymmetry breaking vacua 
which we shall denote \Vi) and IV2). The value of Vi for these minima is given 

by 



1 



3(7ai + 30a2) 



SOMA ± V30^/M'^\^ - 21A/^ai - 90Mfia2 



(6) 



The couplings in this model are chosen in conformity with the requirements 
discussed in [31]. The variables ^ and M in the above equation have dimensions 
of mass whereas A, ai, and 012 are dimensionless. We define the dimensionless 
variable M through the relation M = M^, and choose ai = a2 — 0.1 and 
A = 0.5. The resulting expression for vi becomes 
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(7) 



The value of the dimensionless variable M is varied later in the paper to study 
its effect on stability of vacua containing monopoles. The ± signs in the above 
equation correspond to two distinct vacua, one near the origin due to almost 
exact cancellation and one far from the origin. The vacuum which is far from 
the origin is referred to as iFi). Its properties for AI w 1000 are described below: 



1^1) : Vl 



lOV 



{Vl \V\ Vl) 



loV- 



In contrast, the state IV2) has the following properties: 

m -.vi^ti {V2 \v\ V2) ~ 



(8) 



(9) 



Both these minima have non-zero energy but the energy of IV2) is much smaller 
than that of \Vi). The state IV2) is thus the true vacuum for this model whereas 
\Vi) is metastable. When monopole configurations can exist in the above vacua, 



they shall be denoted \Vi 



monopole\ 



and iVr- 



monopole\ 



respectively. 
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For a given value of V2, the configuration |y™°"°P°'*=^ wliidi asymptotes to 
I^l) at infinity must pass through the field value wi ~ /i corresponding to 1^2) 
near the core of the monopole. This "seeding" of the true vacuum which occurs 
for inside its core can render it unstable as will be shown in section 

15.11 On the other hand, such a seeding does not take place for |y^'"°"°P°''^^ 
and it is thus expected to remain stable. This expectation is confirmed by our 
numerical results which will be discussed in section 15.11 

The phenomenological motivations behind the choice of the superpotential 
given in equation [3] is described in (3l[ . The two terms involving ai and 02 are 
non-renormalizable terms which vanish in the limit M — ?> 00. The desired SUSY 
breaking vacuum of this model turns out to be IV2) and this state survives the 
renormalizability limit. The other vacuum \Vi) gets pushed further away from 
the origin in field space as M increases and eventually becomes tachyonic for 
large enough M. 

3.2. Supersymmetric vacua 

This model contains two supersymmetric vacua when W2 = 0. When this is 
the case, the value of vi for these minima is given by: 



1 



Vl 



2(7ai + 30a2) 



30AfA ± V30v/Af2A2 - 28M/iai - \2QMna2 . (10) 



Once again, the ± signs yield one vacuum close to the origin and one far from 
it. The vacuum further from the origin is denoted as 1^1(5^5^)) for which the 
value of Vl is written as Vi^ . The plus sign in the superscript signifies that 
the plus sign in equation [TU] has been used. Likewise, the vacuum near the 
origin is referred to as |V2(s;7sy)) for which vi is given by . Both these states 
have exactly zero energy, but since Vi < Uj^, a monopole configuration which 
asymptotes into v^ corresponding to \Vi(^susY)) must pass through the field 
value Vl corresponding to \V2{susy)) near the origin. Thus, a "seeding" effect 
similar that of |y™°"°P°''=^ takes place when a monopole is present in \V2(susy)), 
rendering it unstable. This will be discussed in greater detail in section [G] 



4. The monopole ansatz and equations of motion 

We shall set up monopole configurations in the fields Si and S2. The La- 
grangian for the system can be expressed as 

L = -\f;,F^^'' + \{D^^if + \{D^^2? - V^(Si, E2) (11) 

where a = 1, 24 and F^^^ are the gauge field strengths. The covariant deriva- 
tive is expanded as 

7^^E^ = 5^E?-*e[F^,S],]" (12) 
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with i = 1,2 and the adjoints are expressed m terms of the SU{5) generators 
T" as Si — afT°'. We choose a spherically symmetric ansatz for the adjoints 
and the gauge field: 



(13) 



r[l-X(r)], Ag = 0. 



(14) 



To setup the SU{5) monopole, we take the following embedding of SU{2) in 
SU{5): 
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(15) 
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Here, Tq 



1,2,3), and Oa are the Pauli Sigma matrices. The 2x2 
matrices Tq satisfy [Ti,rj] = itijhTk- The SU(5) adjoints Si and E2 given in 
equation [3] are aligned along the generator -^^diag{2,2,2, —3, —3) in isospace. 
In order to obtain the non-trivial topology corresponding to a monopole config- 
uration, we perform a unitary transformation on the fields Si and S2 and also 
the gauge field A to obtain (see [32,] ) 
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(16) 



(17) 



(18) 



where I2 is the 2x2 identity matrix. Using the above results in the lagrangian 
pT|) . the Euler-Lagrange equations of motion become second-order differential 
equations in G{r), H{r) and K{r). Taking e = 1, they are given by 



r'^K" = K{K^ - 1) + K{G'^ + H^) 
d 



r'^G" = 2GK 



r'H" = 2HK' + r'A (l'(G,-f/)) 



(19) 
(20) 
(21) 
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in which the function V{G,H) is obtained by substituting cti = G/r and 
(72 = H/r in the scalar potential ([5]). These equations of motion are ordinary 
differential equations involving only r as a result of our spherically symmetric 
ansatz functions given in (IT51) and (|14p . The function G(r)/r as r ^ co tends 
to a constant value. This means that the function G{r) is proportional to r and 
it therefore diverges at infinity. This is also true for the function H{r). For the 
purpose of our numerical solutions for the monopole, we express equations (jl9l 
- [^T|) in terms of cri and CT2. The result is 

'^'^ \K{K^-l)-K{al+al)^Q (22) 



(Pa2 ^2da2 2 ^ ^2 P^^-Q (24) 
da I r dr r"^ \ 9(72 / 

Equations ([^ - can be solved numerically using relaxation techniques 
after rescaling ai , (72 and the variable r by the VEV vi . A monopole configura- 
tion necessarily implies ai = (72 = at r = 0. From equation (1^^ . the function 
K{r) — 1 as r — > 0. As r — > oo, the function -ftr(r) — ^ while the values of (7i 
and (72 are determined by solving the following set of simultaneous polynomial 
equations: 

-d\r 



(25) 



9(7: 



= 0. (26) 



The initial guess for the solution is chosen to meet the boundary conditions 
discussed above. The value of M is varied to study its effect on the stability 
of the vacua containing monopoles. For a given set of couplings, there are 
two solutions to equations ([^5)) and corresponding to the two translation 
invariant vacua described by and 



5. Stability of the monopole configurations in non-supersymmetric 
vacua 

5.1. Numerical study 

The numerical methods we use are relaxation techniques implemented by 
discretizing the domain over which the solution is required. By dividing the 
interval into a sufficient number of points, we convert the differential equations 
into a set of simultaneous polynomial equations. The initial guess is chosen in 
conformity with the boundary conditions discussed at the end of section 31 

We denote by ^^j^g monopole solution which settles into |Vi) at 

infinity. Likewise, |y^™°"°P°'<=^ denotes the monopole configuration which reaches 
IV2) asymptotically. We have studied the availability of both |y^™°"°P°'^^ a,nd 
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^ymonopo ^^^^ choosing a value of V2 and then varying M. An example of 
^ymonopoie^ when M = 1400 and (T2 = 500 fi is shown in figures |3] and |31 In this 
case, the value of vi is 685. 3/i and this is independent of (72. A similar solution 



exists for IK 



monopole\ 



for which the value of vi is 5.5/i. 




10 15 
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Figure 3: The functions /i(r) = (T-i{r)/v\ and /2('') = o"2('')/fi for ^y^°'^°p°^'^^ with M 
1400 and 1^2 = SOO/j. The value of vi is 685. S/x. 




K(r) . 
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Figure 4: The function X(r) for |y™°"°P°'<=) ^ith ?)2 = 500/i and M = 1400. 
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As discussed in section [31 there is a "seeding" of the true vacuum within 
the core of |y™°"°P°'«^^ jj^ ^j^g above example, the field value ai = vi = 5.5/i 
corresponds to the true vacuum denoted by |V2)- The monopole configuration 
for fji shown in figure |3] must rise from zero and pass through ci = 5.5 fj, before 
reaching its asymptotic value of 685.3/Lt. Furthermore, the energy of the local 
minimum at ai — 685.Sfi increases as M increases. Therefore, for sufficiently 
large values of M, we expect an instability in the configuration 
This expectation is confirmed by the results of our numerical simulations which 
are summarized in figure [5l There is a large region in the parameter space of 
V2 and M for which a numerical solution for cannot be obtained. 



Boundary with monopole ■ 
Boundary without monopoie ■ 



Disallowed without monopole 



Disallowed with monopole 



Allowed 
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Figure 5: The allowed and disallowed regions for The boundary between allowed 

and disallowed regions without a monopole present is defined by equation II27I I. When a 
monopole is is present in \Vi), a large region of the parameter space becomes disallowed since 
^ymonopole^ is no longer available for this region. 



The value of M cannot be increased to arbitrarily high values. There is a lim- 
iting value of M beyond which the metastable vacuum |Vi) becomes tachyonic. 
The condition which must be satisfied for no tachyonic states is 31 1 

2 



V2_ 

Vl 

in which F is defined as 



F 



0L\ + a2 



F ^vf 



— — ai 
M V30 



Vl_ 

AI 



> 



F 



+ C(2 Vl 



(27) 



(28) 



This condition puts an upper bound on possible values of M as shown in figure 
[5l However, when a monopole is present in |Vi), this upper bound becomes 
significantly lowered. 
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The fact that the numerical solution for ^y^°"-°p°'''^'^ cannot be obtained for 
a finite region of the parameter space is indicative of an inherent instability 
in such a configuration. This instability will be studied from a semi-analytic 
approach in the next subsection. The configuration ^^j^ j^g obtained 

for all of the parameter values which were considered. Such a configuration has 
no lower energy state to decay into. 

5.2. Semi- Analytic approach 

The numerical solutions obtained in the previous subsection are all time- 
independent. We can restore time-dependence in the equations of motion and 
study the stability of the solutions without actually solvin g th e time-dependent 



equations of motion. Following the approach discussed in 2Jj , we restore time- 
dependence in equation (f23| : 



When studying the possibility of decaying into jg 

sufficient to study the time-dependence of tii alone. This is because the vacuum 
value of (72 is equal to V2 for both these states. We thus treat a2 as a time- 
independent background field denoted by (T2 (?') in the following discussion. The 
time-dependence of cri(r, t) is decomposed as follows: 

(Ti(r,i) = CTi(r)+p(r)e^"*. (30) 

The function p{r) « ai (r) and ai (r) is the time- independent solution to equa- 
tion (|23p . Substituting equation (|30p in (f29| and linearizing in p, we obtain 



ui'^p - 



^ + P+Pir)]p. (31) 



This equation has the form of a one-dimensional Schrodinger equation with a 
potential 

t/(r)/10V^ = P (MTai^ + 2700^2^ + 1260aia2) (18ct|5-^ + bdf) 
- 20(5Mdf (7ai -I- 30a2) (VsOXai - O.is) 

+ 60I3M^ hx^dl + iX^al - 0.03%/30Acti + 0.0015 



+ mpMdl (7ai + 30^2) [y^QX ai - 0.05 j + — (32) 



where /? — l/450Af^ and ai, (T2, and K are time- independent dimensionless 
functions of r. 

The stability of 1]/"°"°^°'^^ depends on whether or not the frequencies of 
oscillation w in equation pop are imaginary. Real frequencies result in stable 
solutions whereas imaginary frequencies indicate instability. Looking for imag- 
inary modes of oscillation or negative values of uP is equivalent to looking for 
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Figure 6: The equivalent potential for 1^/™°"°?°''=^ with a2 = 500/i. The effect of increasing 
M can be seen to lower the energy of the local minimum near the origin. For large enough 
M, this minimum has a negative energy and the resulting bound state creates an instability 



negative energy bound states of the potential given in equation (|32l) . This po- 
tential is plotted for different values of M with ct2 = 5 for figure 
[6l Notice that there is a local minimum near the origin, and that the energy of 
this minimum reduces as M increases. What this means is that jg 
only stable for small enough values of M. When this is so, there is no negative 
energy bound state possible. When M is large enough, a negative energy bound 
state is possible and the configuration jg longer stable. 



6. Stability of the monopole configurations in supersymmetric vacua 

For the case of the supersymmetric vacua, the monopole configuration set- 
tling into \Vii^suSY)) is denoted \V^su°s'Y{'^ ^ ^^'^ likewise for |^2™sc/sy')'^)- 
turns out that the numerical solutions for \V^su°s'Y{'^ values of M larger 
than 200 can not be obtained. As discussed in the previous section, this failure 
to obtain a numerical solution is indicative of an instability of the configura- 
tion. Thus, even a supersymmetric vacuum containing a monopole can become 
unstable against decay into another supersymmetric vacuum. 

With regards to the stability analyzed from the semi-analytic approach of 
the previous section, the equivalent Schrodinger potential for \V^slJ°s'Y)^) can 
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Figure 7: The equivalent potential for ly^sijgY)^') ■ ^l^^s of M > 200, the negative 

energy bound state of the potential results in an instability in l^/(5'^g^)'^)- 



be derived in the same manner as was done for jg given by 

C/(r)/10V^ = l^t (l47ai + 2700^2 + 1260aia2) 



- 0.01 V30 5-1 A -f 0.0005 



- 4M7CT? (7ai + aOaa) (o.l8 - \/30 Acti) + (33) 



where 7 = 1/90-M^. This potential is plotted in figure[7]and once again, negative 
energy bound states indicate instability. In this case, values of M lying near 
200 represent the boundary between stable and unstable configurations. 



7. Implications for cosmology 

The early universe is characterized by a monotonic reduction in temperature. 
From calculations of thermal effective potentials 33, 34, 3^, we expect a leading 
order temperature correction of the form AT'^af with A > to the effective 
potential of ai . If we absorb this effect in the definition of fj,ef / , we have 



AT'' 



(34) 



If we understand the rescaled quantity M to be expressed in terms of fief f , we 
see that a change in M is equivalent to a change in T. Specifically, we have 
M/fif^ff = M, which shows how an increase in temperature is equivalent to a 
reduction of M. 
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In this way, the increase of M which results in an instabiUty in \y^°'^°P°^'^^ 
can happen naturahy in the early universe with M fixed and the temperature T 
decreasing. Referring to figure [SJ we see that for a fixed value of V2 — ((T2), we 
may start with a high temperature period in the early universe which supports 
the existence of |y^'"°"°P°''=^^ This high temperature is equivalent to low values 
of M . As the temperature drops, the value of M effectively increases and at a 
critical temperature Ti, becomes unstable. Below this temperature, 

the solution \y^°'^°P°^'^^ jg j^q longer available and this region is referred to as 
"disallowed with monopole" in figure [5] As the temperature continues to drop, 
M continues to rise until the state \Vi) becomes tachyonic at a temperature 
T2 < Ti. For temperatures below T2, the translation invariant vacuum \Vi) is 
unavailable and this is referred to as "disallowed without monopole" in figure O 

There is hence an intermediate range of temperatures T2 < T < Ti within 
which \Vi) is available but |y^™°"°P°'^^ jg disallowed. When the universe passes 
through this temperature range, it is possible that causally connected local 
domains settle in the metastable state l^i) instead of the true vacuum IV2). 
Without the presence of monopoles in \ Vi), the transition to IV2) can only take 
place via quantum tunneling. A first-order phase transition of this type can 
result in an inhomogeneous universe 2^, 30j . 

In contrast, if monopoles had formed at higher temperatures in \Vi), the 
resulting states 1^"™°"°?°'*=^ would become unstable within this intermediate 
temperature range. This would result in a classical roll-over to IV2) and prevent 
an otherwise inhomogeneous expansion. In this way, the presence of monopoles 
within \ Vi) ensures a graceful exit from this false vacuum. The universe arising 
out of the type of SUSY breaking model discussed here is hence homogeneous 
like our present day universe. 



8. Concluding Remarks 

The results of this paper share a common theme with the results of [2l|. 
In both metastable vacuum containing a topological defect is shown 

to become unstable against a classical roll-over to the true vacuum. However, 
the exact details of this process depend on the specific model being considered. 



In [2l|, the model of gauge mediated supersymmetry breaking which we had 
studied permitted the existence of cosmic strings. These string configurations 
were responsible for the seeding of true vacuum bubbles within their cores. 
Furthermore, the metastable vacua which became unstable in the presence of 
the strings were the phenomenologically desired SUSY breaking minima. The 
true vacuum was an undesirable minimum in which SU{3)c was broken. In this 
way, the cosmic strings played a potentially disastrous role in the model. 

The role of the topological soliton is quite the opposite and much more 
benign in the present case. Here, the SU{b) grand unified group supports the 
existence of monopoles and as we have discussed, the metastable vacua are 
phenomenologically undesirable. There is hence the need for a graceful exit 
from the false vacuum and this is where the monopoles come in. If the universe 
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is trapped in a false vacuum and decays by quantum tunneling, the resulting 
universe can be inhomogeneous. The type of classical monopole instabilities 
discussed here ensure that the metastable vacuum can become unstable and 
undergo a "roll-over" to the global minimum. This effect hence provides a 
graceful exit from a potentially inhomogeneous expansion. 

In general, any model based on a gauge group which supports the existence 
of topological solitons can undergo the process described here. Whether or 
not the metastable vacua containing solitons are desirable depends on the given 
model. Hence, the stability of a metastable vacuum against quantum mechanical 
tunneling is not always a sufhcient condition to ensure the viability of a model 
containing solitons. 
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